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ABSTRACT 

Context. We analyze 134 ks Chandra ACIS-I observations of the Galactic Centre (GC) performed in July 2011. The X-ray image with 
the field of view 17' x 17' contains the hot plasma surrounding the Sgr A*. The obtained surface brightness map allow us to fit Bondi 
hot accretion flow to the innermost hot plasma around the GC. 

Aims. Contrary to the XMM-Newton data where strong 6.4 keV iron line was observed and interpreted as a reflection from molecular 
clouds, we search here for the diffuse X-ray emission with prominent 6.69 keV iron line. The surface brightness profile found by us 
allowed to determine the stagnation radius of the flow around Sgr A*. 

Methods. We have fitted spectra from region up to 5" from Sgr A* using a thermal bremsstrahlung model and four Gaussian profiles 
responsible for emission lines of Fe, S, Ar, and Ca. The X-ray surface brightness profile up to 3" from Sgr A* found in our data 
image, was successfully fitted with the dynamical model of Bondi spherical accretion. 

Results. We show that the temperature of the hot plasma derived from our spectral fitting is of the order of 2.2-2.7 keV depending 
on the choice of background. By modelling the surface brightness profile, we derived the temperature and number density profiles in 
the vicinity of the black hole. The best fitted model of spherical Bondi accretion shows that this type of flow works only up to 3" and 
implies outer plasma density and temperature to be: = 18.3 ±0.1 cm“^ and 7°“* = 3.5 ± 0.3 keV respectively. 

Conclusions. We show that the Bondi flow can reproduce observed surface brightness profile up to 3" from Sgr A* in the Galactic 
Center. This result strongly suggests the position of stagnation radius in the complicated dynamics around GC. The temperature at 
the outer radius of the flow is higher by 1 keV, than this found by our spectral fitting of thermal plasma within the circle of 5". The 
Faraday rotation computed from our model towards the pulsar PSR J1745-2900 near the GC agrees with the observed one, recently 
reported. We speculate here, that the emission lines observed in spectra up to 5" can be interpreted as the reflection of the radiation 
from two-phase regions occurring at distances between 3 - 5" of the flow. The hot plasma in Sgr A* illuminated in the past by strong 
radiation field may be a seed for thermal instabilities and eventual strong dumpiness. 
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1. Introduction 

The central region of the Milky Way hosts a supermassive black 
hole (BH) with we ll established mass Mbh = 4.4 x IO^Mq 
( iGenzel et al.ll201(ill . at a distance R - 8.4 kpc. The surround¬ 
ings of the black hole are uni que. It contains a number of rel¬ 
atively cold molecular clouds (IZhao et al.ll2009ll located within 
15' from Sgr A*, and a hot p lasma of temperature around 2 keV 
situated within a few arcmin (iBaganoff et al.ll200^ . Both phases 
of gas show X-ray emission with prominent iron Kq, line, which 
may arise due to reflection of hard radiation emitted at the very 
center of the Galaxy nucleus. 

The cold phase emits neutral Fe Ka line at ener gy 6.4 keV, 
as in c ase of molecular clouds reported recently bv iPonti et alJ 
(I 2 OIOI) . These cold clouds situated at the distance of 30-60 pc 
from the Galactic Centre (GC), are well resolved by the XMM- 
Newton satellite, enabling proper studies of their velocities and 
the nature of the reflected continuum. Tracing history of fluores¬ 
cent reflection, it is possible to c apture historic luminosity vari¬ 
ations of Sgr A* dYu et al.ll201lll . 
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The hot phase emits the Fe Kq, line at energy 6.690 keV, in¬ 
dicating that iron atoms are in helium-like ionization state. X - 
ray emission from a diffu se hot gas was observed in extende d 
regions around the GC dPark et al.l l2004t iMuno et ^ l2004ll . 
Particularly, it appears to be concentrated within the central 2-3 
pc eastwards from Sgr A*, named Sgr A East, and most proba¬ 
bly is caused by a young supernova remnant (SNR) dMaeda et al.l 
I 2 OO 2 I) . Additionally, the hot phase lies on the line of sight to¬ 
ward Sgr A* and it may accrete on supermassive black hole in 
the form of spherical Bondi accretion or radiatively inefficient 
accretion flow (RIAF). The dynamical model of the accretion on 
the GC is still not well established since it is clearly seen that this 
region may be a mixture of hot p lasma with cold clouds i.e. mini¬ 
spiral region dZhao et al.l 1201 Oh . The accretion of hot plasma 
fe d by stars due to mass and ene rgy exchange was proposed 
bv IShcherbakov & Bagano^ d2010h to expl ain the observed sur- 
face brightness profile. On the other hand, ICzernv et al.l d2013h 
constructed a model of multiple accretion events caused by cold 
clouds to trigger the Sgr A* activity. Such hot and cold two- 
phase medium can be formed due to thermal instability caused 
by irradiation with hard X -rays from the GC, as recently shown 
bv iRozariska et alJ d2014ll . Thermally unstable clouds can be lo- 
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cated close to the center at 0.008-0.2 pc, i.e. within 5", and they 
can survive hundreds years. The radiation reflected from those 
clouds can produce hot iron line in that region. 

Since the extension of the hot gas is an order of magni¬ 
tude lower than of the cold phase, Chandra and Suz aku satellites 
are better for detailed studies of those regions ( e.g. iMaeda et alJ 
I 2 OO 2 L iKovama et al.l 120071: IWang et all 20131 and references 
thereinl. Recently. IWang et al.l (120131) have reported results of 
spectral htting of the longest 3 Msec Chandra observations of 
the 1.5" circular region around Sgr A*, and the spectrum pre¬ 
sented in this paper is the most detailed ever published. But the 
surface brightness prohle from dynamical model is not htted in 
this paper. The authors ht model spectrum computed as a ther¬ 
mal emission from the distribution of temperatures and densities 
according to the analytical formula appropriate for RIAF (see 
Sec. l4.2l in this paper. 

The purpose of this paper is to study the hot plasma around 
Sgr A*. Analysing X-ray image we aim to put constraints on 
the physical parameters of the hot phase, which may be an im¬ 
portant material to be a source of radiation from Sgr A*. We 
present archival Chandra ACIS-I observations of the GC made 
in July 2011 with the field of view 17' x 17' containing both the 
Sgr A* and Sgr A East. The data were collected with three time 
intervals, very short one after another. Therefore, spectral and 
imaging coadding was very accurate. 

From Chandra observations, we extracted spectrum from 
one circular region around Sgr A* with radius of 5". We per¬ 
formed detailed spectral analysis by htting a thermal plasma 
with a strong K(, line from helium-like iron. Additionally, 
we htted other emission lines seen in residuals form C, Ar, 
Ca, but the detection o f last two was marginal. Following 
IShcherbakov & Baganohl (l2Q10li . we extracted surface bright¬ 
ness prohle around Sgr A* within the radius of 10" with sub¬ 
pixel accuracy. We htted the dynamical model of the hot Bondi 
accretion how by iMoscibrodzka et ^ (l2009l) to this imaged 
emissivity prohle. To make a ht, for each computed model we 
constructed the map of emission around GC convolved with 
Chandra exposure map of ACIS-1 chip. In our observations 
Sgr A* is located off-axies, and we constructed point spread 
function (PSF) to account for this particular observation. Final 
models were convolved with normalized PSF, before comparing 
to the data. We point out here, that only basic spectral htting is 
done in this paper due to the short time exposure. In contrast to 
IWang etalj|2()13li . we do not compute the line emissivity from 
our model. Instead, we foc us on modeling the sub-pixel su rface 
brightness prohle (see also lShcherbakov & Baganoftll2010h . 


The main result of the observed surface brightness prohle 
modelling is the derivation of electron density and temperature 
prohles of the How from 3" down to the black hole horizon. The 
best htting model is for temperature and density at the outer How 
radius T™' = 3.5 + 0.3keV, andn°“' = 18.3 + 0.1 cm“^. Ourhtis 
not valid outside 3" indicating the location of the stagnation ra¬ 
dius of the dynamical how around GC. Ouside stagnation radius, 
m atter can po s sess o uthow as it was indicated in RIAF model 
of IWang et al.l (l2013h . The temperature at the outer how radius 
is in good agreement with that derived from spat ially resolved 
spectral modelling presented in I Wang et al. ( 2013h. and then re¬ 
jected by RIAF model. In addition, Wang et alJ (l2013h consid¬ 
ered RIAF model, where the line emissivity was computed from 
the predicted RIAF temperature prohle. Such line htted to the 
data indicated the outer plasma temperature to be 1 keV, which 
is lower than in our model. In this paper we compare temperature 
and density prohles implied for both models, RIAF and Bondi. 


Table 1. Chandra satellite observation log. 


ObsID 

Date 

Exp. time 
[ks] 

GTI 

[ks] 

12949 

2011-07-21 

58.48 

48.88 

13438 

2011-07-29 

66.18 

59.79 

13508 

2011-07-19 

31.49 

26.10 


Fitting the canonical model of thermal bremsstrahlung plus 
Gaussian prohles for all lines indicates temperatures for the 
plasma around Sgr A* in the 2.2-2.7keV range, depending on 
the choice of ba ckground. And agai n, those values are lower than 
that obtained bv lWang et al.l (120131) with the same model for the 
1.5" Sgr A* spectral htting, but higher then 1 keV temperature 
derived by them using RIAF emissivity prohle. The EW of the 
iron line we hnd is in the 0.9 - 1.2keV range. Within the uncer- 
tai nties, our EW of th e iron line is consistent with that estimated 
by IWang et al.l (1201 3h at the Icr conhdence level. Since the re¬ 
sult is degenerated according to plasma temperature and warm 
absorption, the broad-band spectra are needed to discriminate 
between models. 

In addition, we calculated the Faraday rotation (RM) to¬ 
wards the Galactic Center and towards the pulsar located at 3" 
in the projected distance a way from the GC. T he hrst result 
matches with observations (iMarrone et aLlKOO^ only with ad¬ 
ditional assumption of the very weak magnetic held, while the 
second result agrees with the recently observed value very well 
(lEatough et al.ll201^ . 

Our paper shows that Bondi accretion can well represent the 
hot phase accretion on the GC up to 3". We argue that the radi¬ 
ation originating from the very close vicinity of Sgr A*, illumi¬ 
nating external zones of the Bondi how may produce two-phase 
medium when all lines are created. Th e exact distance of suc h 
clumps and their lifetime are given in iRozariska et alJ (l2014l) . 
The temperature of such clouds varies from log(r) = 5.5 - 7.5 
in Kelvins, and the warm clouds have longer lifetimes than the 
cold clouds. We postulate here, that they can be responsible for 
the formation of iron line at energy 6.7 keV, seen in the Chandra 
data. Better data are needed to fully conhrm this hypothesis. 

The structure of the paper is as follows; in Sec.|2]we present 
the details of data extraction that resulted in obtaining images, 
brightness prohle and spectra of hot plasma in the neighbour¬ 
hood of the GC. Sec.[3]contains results of spectral htting, while 
Sec. |4] shows results of surface brightness prohle htting in the 
vicinity of Sgr A*. Discussion and conclusions are presented in 
Sec.|5] 


2. Observations 

The Chandra X-Ray Observatory has observed the central re¬ 
gions of the Galaxy multiple times. For our analysis, we have 
selected three observations carried out in July 2011. the longest 
unpublished observations available in the Chandra archive. They 
were performed for a very similar period of time, allowing to 
avoid signihcant changes with time and astrometry. Table [1] 
presents identihcation numbers, starting dates and exposure 
times of the observations. All of them were done using ACIS- 
I spectrometer with the held of view 17’ x 17’. It corresponds to 
the region of 40 x 40 pc in the center of our Galaxy. The point¬ 
ing coordinates were a = 17^45^59.^70, <5 = -28°58'15."9, 
and were 4.'8 away from Sgr A*. 
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to minimize the influence of the molecular clouds variability we 
i. 5 oe-o 6 made two tests: i) point sources were removed from the ACIS-I2 
CCD chip, where diffuse emission was marginal, filtering counts 
iiirts to the 0.5-8 keV energy band ii) photons with energy from 7 
to 10 keV in the whole ACIS-I field of view were considered. 
Biie -07 In both cases similar good time intervals (GTl) were obtained. 

Hereinafter, we decided to use GTl found in case (i), as they set 
594 e.o 7 more stringent limits. For all observations the light curves are of 
good quality and the background flaring is insignificant. Finally, 
4.3oeo7 we choose GTI= 134.77 ks, which covers 86% of the total expo¬ 
sure time. 

3 0494)7 Although during all observations the telescope was pointed 

at the same target, we have found a small shift between three im- 
2.07M7 ages (a few arcsec). In order to correct the absolute astrometry 
of the data, we searched for point sources using the wavdetect 
1 2894)7 tool. 53 bright sources were found and their coordinates were ob¬ 
tained from the SIMBAD astronomical database^. Subsequently 
6 13908 vve have applied minor corrections to the WCS using the repro- 
JECT_ASPECT script, which calculates shift between the wavdetect 
output and the absolute coordinates from a catalog. 
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Fig. 1. Flux image of the GC in the 0.5-8 keV energy range (up¬ 
per panel). Sgr A* is ma rked with a black cro ss. A few bright¬ 
est filaments reported by Llohnson et alJ (l2009l) are marked with 
white (small) ellipses and named following the names given in 
previous detection . The green (large) ellipse shows the region 
used for background spectrum extraction (the case of distant 
background, see Sec. 12.41 for details). Bottom panel shows flux 
image of the GC in the energy range of 6.3-7.0 keV, where iron 
line is created. Emission presented here is the Kq, line emission 
after subtracting the continuum flux. The flux is presented in as- 
inh scale, since it behaves well for negative count rate values 
which might originate during subtraction process. 

2.1. Data reduction 

All subsequent analysis was performed using CIAO 4.4 soft- 
war41 and CALDB version 4.5.3. 

We started our analysis by processing individual observa¬ 
tions. The corrections for charge transfer inefficiency, bad pixel 
removal, and background flaring were done in the whole energy 
band, 0.5-8.0 keV. The lc_sigma_clip tool was used to clip the 
data that deviates from the mean by more than 2 sigma. In order 

' http://cxc.harvard.edu/ciao/index.html 


2.2. Images in continuum and in iron iine bandpass 

We have regraded event files to a common tangent point with 
REPRorECT_EVENTS and we have merged them using dmmerge. 
Similar procedure was applied to the images, beside that the re- 
projectjmage and dmimgcalc scripts were used instead. We have 
created exposure-corrected image in the 0.5-8 keV energy band 
using the fluximage script. The flux image is presented in Fig.[T] 
upper panel, with the position of Sgr A* marked by the black 
cross. At this figure we h ave also indicated and named the bright¬ 
est filaments reported bv iJohnson et al.l ( l2009h . Nevertheless, our 
total exposure time was too short to make spectral analysis of 
these filaments or discover new ones, so we do not proceed with 
their analysis. However, by merging the three observations we 
obtained enough counts to study the hot extended plasma around 
Sgr A*. Since Chandra X-ray observatory has the best spatial 
resolution among all currently working satellites, in Sec.|4]we 
present analysis of the surface brightness profile around Sgr A* 
with 1" scale resolution. 

Since we were interested in the spatial distribution of the 
ionized iron Kq, line, additional images were created in the 4.5- 
6.3 keV (the continuum) and 6.3-7.0 keV (broad line) energy 
bands. The former was subtracted from the line image with a 
normalization factor of 0.12. This constant was calculated on the 
assumption that the continuum emission is modelled by an ab¬ 
sorbed power-law function with photon index F = 3 and hydro¬ 
gen column density Nh - 1 0^^ cm~^. The contin uum spectrum 
was modeled using Sherp^l (iFreeman et al.ll2()()il) and simulated 
with the FAKE_PHA command. The response matrix and effec¬ 
tive area files were created with mkacisrmf and mkwarf (sample 
RMFs and ARFs were also taken from CALDB, similar value 
was obtained). The normalization factor is simply the ratio of 
the flux in the line to the flux in the continuum energy band. 

The image in iron line emission after subtracting continuum 
image is presented in the bottom panel of Fig. [1] The majority of 
point sources and diffuse emission around Sgr A* are properly 
subtracted. Some faint molecular clouds remained after subtrac¬ 
tion because the e mission from the cold Fe line overlaps broad 
line energy band (iPonti et al.l 1201 Oh . We have constructed the 
flux image in line, but the emission was negligible. 

^ http://simbad.u-strasbg.fr/ 

^ http://cxc.harvard.edu/sherpa4.4/ 
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Fig. 2. X-ray image of the 10" region around Sgr A* (circle). 
Surface brightness profile was extracted integrating counts over 
annuli from 0 to 10". The influence of point sources was rejected 
as shown in the image. 



Fig. [T]bottom panel is presented in the asinh (the inverse hy¬ 
perbolic sine) scale. In case of the high count rate this scale is 
practically logarithmic; for low values it changes linearly. It is 
especially useful in subtracted images i.e. the iron line, as it be¬ 
haves well for negative count rate values which might be present. 

2.3. Sgr A* surface brightness profile 

The surface brightness profile can be constructed in counts 
per pixel squared as a function of distance from the 
BH. It was previously done for much longer observations 
(IShcherbakov & Baganoffll2010h with the exposure time 953 ks. 

Here, we repeat this analysis for our data, but we fit them 
with a different model (see below Sec. |4|. The size of Chandra 
pixel is 0.5", but the position of the satellite over the duration 
of observation can be determined with the 0.1" accuracy by 
comparin g with the known position s of bri ght point sources. As 
noted by Shcherbakov & Baganofil ( 1201 Ol) we can achieve the 
pixel resolution accuracy in the surface brightness profile from 
knowing the orientation of the detector pixels at the given time. 

In Fig. |2] we show pixel image of the close vicinity of the 
GC. It was necessary to exclude some point-like features in the 
Sgr A* neighbourhood because they broke spherical symmetry 
around the source. In order to obtain a radial profile, the net 
counts in a set of concentric annuli around Sgr A* were mea¬ 
sured with DMEXTRACT tool. The width of annuli is 0.25", the 
radius of the outer annulus is 10", and is marked in Fig.|2]by a 
black circle. 

The radial profile is monotonically decreasing in the whole 
range of radii. The final profile of the surface brightness with 
errors is presented in Fig.[2by red points. 

To make any further analysis of the brightness profile on the 
sub-pixel scale, we have to account for the point spread func¬ 
tion (PSF) appropriate for the Chandra X-ray telescope. Any 
model fitted to the profile has to be convolved with the PSF. 
To construct the ACIS-I PSF for our observation we used the 
Chandra ray-tracing program, ChaRlQ, to simulate the Sgr A* 
(point source) photons scattered by the Chandra mirrors, and the 
MAR}^ software version 4.5.0 to projects the simulated rays 

http://cxc.harvard.edu/chart/runchart.html 
^ http://cxc.harvard.edu/chart/threads/marx/ 



Distance from BH, arcsec 


Fig. 3. Vicinity of the black hole up to 10" from Sgr A*. Red 
points with errors represent observed profile, and are connected 
by continuous line. Black long dashed line is normalized PSF. 
Two models convolved with PSF are also presented: green dotted 
line for T™' = 1 keV, while blue short dashed line for = 3.5 
keV. Residuals are shown for best fitted model in the range up to 
10", i.e. for T™' = 3.5 keV. 


onto the detector plane. The MARX output was used to extract 
the PSF profile for convolution with our model spectra. ChaRT 
requires the position of the point source on the chip, exposure 
time and the point source spectrum to run the simulations. For 
these we used the position of Sgr A* in our stacked observa¬ 
tion, total GTI exposure, and the best fitting model to the inner¬ 
most region of Sgr A* (5”, Table |2]l. The normalized PSF con¬ 
structed for our particular observation is presented in Fig. [3]by 
black cr osses and long dashed line. The P SF is wider than in the 
paper bv IShcherbakov & BaganolS (1201 Oh . due to the fact that in 
our observations Sgr A* was slightly off-axis. 

2.4. Extraction of the spectrum 

We have extracted spectra from regions where the emission 
of the line was the highest, especially from the region around 
Sgr A*. Sgr A East region is treated separately and it will be 
presented in the forthcoming paper. 

For the purpose of this work, spectrum from circular region 
of the radius of 5" around Sgr A* was extracted and analysed 
in details. This was done to see how far from the BH the iron 
line emission is crucial, and to estimate the temperature at the 
outer radius of the dynamical flow, which is independently de¬ 
rived while the surface brightness profile is fitted as presented in 
Sec.l4Tl 

Prior to creating the spectrum, a few point-like sources were 
excluded since they might contaminate the results. The spectrum 
was obtained with the specextract script for each observation 
independently and merged with the combine_spectra tool. We 
consider two different regions for background: a distant region 
marked as a green ellipse in the upper panel of Fig. [1] and a local 
region defined as an annulus between 6" - 18" around the 5" 
source region. The distant background region does not contain 
any point sources, and the diffuse emission is negligible inside 
the ellipse. In the case of local background the possible point 
source contamination was investigated and corrected for during 
the data reduction process. After background subtraction, the 5" 
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Table 2. Best fit parameters for the 5" region around SgrA*. We present two cases modelled with two different statistics: CSTAT 
for the case where source and background are fitted simultaneously, and;^^^ for background subtracted analysis. In the second case 
we consider two background regions: a distant background marked with the green ellipse in the upper panel of Fig.[T] and a local 
background defined as an annulus between 6" - 18" around the 5" source region. In the first row we show results obtained with 
a model consisting of thermal bremsstrahlung with Galactic absorption (tbabs) and Fe Gaussian line only. The second row shows 
model parameters when three more Gaussian lines are added to account for S, Ar, and Ca emission. Both fits were done using CSTAT 
statistic. The third and fourth rows describe background subtracted analysis (using the local and distant backgrounds, respectively) 
performed on data grouped according to S/N = 3 using the;^^^ statistic. 


Bkgr. 

Nh 

kT [keV] 

Line 

Ek„ [keV] 

Line 

Ek. [keV] 

Line 

Ek„ [keV] 

Line 

Ek„ [keV] 

Stat. 

d.o.f. 

Reg. 




EW [keV] 


EW [keV] 


EW [keV] 


EW [keV] 

Val. 


Distant 

7 

'•^^-0.15 

Q ^7+0.08 

-0.09 

Fe 

6-728-:-= 


— 


— 


— 

CSTAT 






0.58 ± 0.06 







1759 

1019 

Distant 

^•^‘’’-0.25 

^•^^-0.07 

Fe 

f\ 79Q+0.064 

''^^-0.042 

S 

2.480!--- 

Ar 

3 153+0-05 

Ca 

3.861 

CSTAT 






^•^^-0.24 


0.18!«:“ 


0 05 
^•^‘-^-0.02 


^•^“^-0.01 

1254 

1008 

Local 

^^•-’^-0.39 

-0.07 

Fe 

6 745+U.118 

VJ. /‘t^-0 069 

S 

z,.‘t07_oo3o 

Ar 

3.215::: 

Ca 

3.865::: 







1 1Q+0.89 
^•^^-0.48 


0.20!“:°’ 


u.uU-003 


0 05 

^•^‘-^-0.03 

117.32 

208 

Distant 

Q 4^+0.31 

9 ^Q+O.U/ 

^•-’0-0 09 

Fe 

f\ 79 1 +0.004 

'■^4_o_o45 

S 

2.488!-- 

Ar 

3 140+0-^59 

j.i^U-0 040 

Ca 

9 86^/^+0 055 

■J-ot)O_0,086 







0 92+0-32 
-0.26 


0-21!^ 


^•^^-0.03 


0 04+0 03 
-0.02 

122.64 

198 


source photon’s number in the 0.5-8 keV band is equaled to 2937 
and 2358, for case distant and local backgrounds, respectively. 
Full description of our spectral analysis is presented below in 
Sec.|3] 

3. Spectral fitting 

The spectral fitting of the data was performed with the Sherpa 
4.4 fitting package. In all plots presented in this section, the data 
were grouped requiring the signal-to-noise ratio S/N - 5, for 
presentation purposes only. All errors on the spectral parameters 
correspond to the 68% confidence level (Icr) for one significant 
parameter. 

In the first step, considering case A (distant background), 
we rejected simple background subtraction process. Instead, the 
source region and corresponding background were modelled si¬ 
multaneously using the CSTAT statistic. We used a phenomeno¬ 
logical model composed of a thermal bremsstrahlung emission 
(bremss) to describe the continuum, and a Gaussian profile to 
account for the iron Kq, line emission. The total model was mul¬ 
tiplied by warm Galactic absorption (tbabs). The background 
data were parametrized with a power law model. Resulting fit 
parameters are presented in Tab.|2]in the first row. Figure |4] (up¬ 
per panel) shows the data and the best fit model together with 
residuals. The fit resulted in CSTAT=1759 for 1019 degrees of 
freedom (d.o.f.). The residuals below 4keV hint at the presenc e 
of additional emission lines detected e.g. bv lWang et al.l (l2013h . 

In the second step, using the same source/background re¬ 
gions statistics, we add three fit Gaussian profiles to t he model to 
accoun t for S, Ar, and Ca emission lines reported bv lWang et alJ 
(l2Q13ll . The widths of these additional lines were fixed at cr = 
10“^ keV. Certain residuals around the energy of the S line were 
visible also in the background data, and thus the S line was added 
also to the background model (again with cr fixed at 10“^ keV). 


These changes resulted in CSTAT = 1696 for 1007 d.o.f., a de¬ 
crease in the temperature of the plasma, and increase in the col¬ 
umn density of the absorber. 

Further improvement to the fit quality was obtained by ad¬ 
justing the curvature of the background model, i.e. replacing 
the simple power law model with a 3rd degree polynomial, and 
adding a Gaussian with negative normalization to account for a 
feature visible below 1 keV. These updates did not affect the val¬ 
ues of plasma temperature nor the column density, but they led 
to a substantial decrease in the fit statistics, CSTAT = 1254 for 
1011 d.o.f. The best fit parameters are shown in the second row 
of Tab. 12 The data and best fit model are presented in Fig |4] 
lower panel. 

We investigated if the change in the values of the Ah (from 
7.44;^Q[g X lO^^cm^^ to 9.15);g25 x lO^^cm^^) and kT (from 
3.57);QQgkeV to 2.66keV) could be due to a degeneracy 
between these two parameters. Indeed we found that at each 
step we are able to find both a 1ow-Ah high-kT solution, and a 
high-An low-kT solution with the statistics that differed only by 
ACSTAT = 5-32 (for d.o.f. changing from 1019 to 1007, corre¬ 
spondingly), with the high-An low-kT solution always with the 
lower CSTAT. We interpret it as an indication that these two pa¬ 
rameters are indeed degenerate in our modelling, with a slight 
hint towards the high-An low-kT solution. 

Tab. |2] shows also the the best-fit model parameters for two 
additional fits with total model tbabs -i- bremss -h 4 lines, us¬ 
ing the default Sherpa statistic (chi2gehrels), and background 
subtracted data grouped requiring that S/N — 3. We performed 
these fits in order to check how different ways of treating the 
background affect our analysis. The third and fourth rows of the 
table presents results obtained using the local and distant back¬ 
grounds, respectively. The corresponding plots are presented in 

Fig. 13 
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Fig. 4. Modelling of the spectrum extracted from the 5" cir¬ 
cular region around the Sgr A*. Both panels present hts with 
distant background. The upper panel shows a model consisting 
of Galactic absorption (tbabs), thermal bremsstrahlung and one 
Gaussian line to account for the iron emission. The lower panel 
shows a model with three additional Gaussian lines for the S, Ar 
and Ca emission. 


The change of the statistics does not change the results if 
the same (distant) background is used. All model parameters are 
consistent within their uncertainties. In the case with the local 
background, we observe a further decrease in the plasma tem¬ 
perature to a value kT - 2.24;|;qq!^ keV and increase in the ab¬ 
sorbing column. Ah = 10.52);Q3g x lO^^cm^^; these values are 
not consistent at 2cr level with the estimates from the other mod¬ 
els including all 4 emission lines. 

We conclude that the continuum is well fitted by a thermal 
bremsstrahlung model indicating temperatures for the plasma 
around Sgr A* in the 2.2-2.7 keV range, depending on the 
ch oice of background . These values are lower than that obtained 
by I Wang et aljl2013h with the same model for the 1.5" Sgr A* 
spectrum; they quote 3.5(3.0,4.0)keV. The EW of the iron line 
we find is in the 0.9- 1.2keV range, but its uncertainties are high, 
due to the uncertainty in deriving the relative contribution of the 
continuum and line, the number of model parameters, and the 
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Fig. 5. Results of spectral htting of 5" region - upper panel for 
distant background, lower panel for annulus 6" - 18". The total 
model marked with a solid red line is a sum of Galactic absorp¬ 
tion (tbabs), thermal bremsstrahlung and four Gaussian lines for 
Fe, S, Ar and Ca emission. 


quality of the data. Wit hin the uncer t aintie s, ourEW is consistent 
with that estimated by IWang et al.l (1201 3h at the Icr conhdence 
level (EW= 691(584,846) eVs). Our data are of not sufficient 
quality to study the Fe line properties (e.g. its evolution with the 
distance from Sgr A*) in details and for this purpose we address 
the rea der to the wo r k oflWang et alJ (l2013h . The absorbing col¬ 
umn in I Wang et al.l (1201 3i) . 10.1 (94.11.11 x lO^^ cm is also 
consistent with our results at the Icr conhdence level. Note, how¬ 
ever, that this value of Ah is inconsistent with the absorbing col¬ 
umn they hnd for the Haring Sgr A* spectrum, and based on this 
they tend to discard this phenomenological Htting as unphysical 
and continue with the development of their RIAF model, where 
they Hnd plasma temperature ~ 1 ke V. Our value of kT is some¬ 
what intermediate between those of lWang et al.l (1201 3h resulting 
form their phenomenological and RIA F models. Howeve r, in the 
following, we show alternative to the IWang et ^ (1201 3h Htting 
of the surface brightness proHle. Instead of computing emissiv- 
ity from the gas with temperature and density proHle analytical 
expressions , whic h mimic RIAF dynamical model as done by 
IWang et al.l (1201 3h . we Ht dynamical Bondi flow directly to the 
luminosity map, without any proHle functions. 
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4. Hot flow around Sgr A* 

Chandra ACIS-I camera spatially resolves close vicinity of 
the central supermassive black hole and the observations con- 
strain the model of hot plasm a around Sgr A*. Previously, 
IShcherbakov & BaganoffI (1201 Ob analyzed 953 ksec Chandra 
data and they have htted model of hot plasma interacting with 
stars up to 5" from the Sgr A* complex. They have provided 
density and temperature radial prohles of the hot plasma at the 
outer radius of the fl ow to be equal: T - 2 keV and We = 60 
cm“^. More recently. IWang et ^(l2013h reported 3 Msec obser¬ 
vations of Sgr A*, where they have concluded that RIAF model, 
when the spectrum is htted by calculating emissivity from RIAF 
temperature and density prohles, provides an excellent spectral 
ht indicating the plasma temperature to be 1 keV. Unfortunately, 
the authors did not ht surface brightness prohle to the data. 

In the data extracted by us, the collected counts per pixel 
square are above 30 in the center and a few at 10" (Fig. |3 red 
line). Below we present the ht to the surface brightness prohle 
with the model of the hot Bondi how. The purpose of our studies 
is to show to which extend the Bondi how works. In principle 
in the very centre star inhuence can be inefficient as discussed 
by IShcherbakov & BaganoffI (l2010h . In any of such models or 
RIAF model with outhow, there exists a stagnation radius, which 
indicates the border between the matter inhowing or outhowing 
from the very center. At the stagnation radius the Bondi how 
stops to be valid, as we show in the section below. 


4. 1. Dynamical model of the hot accretion flow 

We compute model brightne ss prohle bas ed on a standard spher¬ 
ical accretion how theory dBondil [19521) . The plasma density, 
temperature, and velocity radial prohles are obtained by solving 
the general Bondi equations (the conservation equations given in 
IShapiro & Teukolskvl 19831 equations G.21 and G.22). In partic¬ 
ular we solve the following Bernoulli equation: 


h 





1 - 3 * cs? 


( 1 ) 


where Huid enthalpy is h = 1/(1- cs^Kjad - 1)), cs is the speed 
of sound, u’^ is the huid velocity, and the right-hand-side is the 
integration constant, found by taking the h and cs at the sonic 
radius, denoted with underscore s. The speed of sound at the 
sonic radius r, is given by C5^ = 1/(2 * rj - 3) (see Eq.Gl? in 
IShapiro & Teukolskvlll98.3h . The adiabatic index is assumed to 
be yad ~ 5/3. Although we use equations that generalize Bondi 
equations for the Schwarzschild spacetime, our solutions at large 
distances from the black hole naturally converge to the classical 
Bondi solution. 

The above Bondi model has two (except yad which is hxed 
in our model) free parameters: density and temperature of 
plasma r°“‘, both at the outer radius Rout ^ 10". These vari¬ 
ables are assumed far away from the black hole, and so they are 
equivalent to the “inhnity” values (i.e. n™' s noo, T™' = Too). 
The mass accretion rate onto the black hole in the Bondi models 
in the units of g s ' is: 


M = TttT, 


(GMbh)Voo 

*^S,co 


( 2 ) 


where Ts is a function of adiabatic index yad (for our yad ~ 5/3, 
ds = 0.25). We integrate the model equations from Rout down to 
the black hole horizon. The list of models explored in this work 
is given in the Tab. [3] 


To create synthetic brightness profiles, we first generate im¬ 
ages of the inner region of the GC 20" x 20" using a ray trac¬ 
ing method. Our scheme performs radiation transfer through the 
three dimensional spherically symmetric accretion flow model 
and records the radiation flux density at each of the detector pix¬ 
els. In our modelling, we assume that the X-ray radiation is pro¬ 
duced by bremsstrahlung. We solve the radiation transfer equa¬ 
tion along many rays, using the following bremsstrahlung emis¬ 
sivity function: 

efF,y = 6.8 X 10^^^ nl exp j. (3) 

The plasma absorptivity is also taken into account and it is 
derived from the Kirchhoff’s law of thermal radiation: ag^y = 
eg^ylBy, where By is the Planck function. 

To convert the theoretical flux to counts detected by ACIS-I 
Chandra CCD, the model fluxes are convolved with the instru¬ 
ment response. The total number of counts detected during the 
exposure time Tgxp is: 

Ctot = J'rexpR’vA effiv)v ^ exp{-abs(v)NH)dv, (4) 

where visa frequency corresponding to energies from 0.5 up to 8 
keV. Additionally, we assume that absorption also takes place in 
the ISM between the outer boundary of accretion model and the 
detector. Hence, we multiply the flux by additional exponential 
function, and Ah = 10 ^^ cm“^ is a typical column density of 
the ISM. Effective area of the instrument Aeff(v) and absorption 
coefficient in the ISM abs (v) are energy dependent; we ado pt the 
same dependencies as in IShcherbakov & BaganoffI (l2010l priv. 
communication). 

The final modelled emissivity profiles are obtained by map¬ 
ping the square images (20" x 20") onto 40 concentric rings (the 
same rings that were used for extracting the data. The total count 
rate collected on each ring is divided by that ring area (ring area 
is given in a units of pixels). 

Thus computed profiles have to be convolved with the PSE 
derived for our particular observation in Sec. 12.31 Examples of 
convolved models of Bondi hot accretion are shown in Eig.|6]by 
dotted and short dashed lines for two outer temperatures equal 
1 and 3.5 keV respectively. Eor fitting procedure, we have com¬ 
puted a grid of models with different temperatures and densities 
ranging in T™' = 1-16 keV and =1-60 cm“^ (see Tab.[3]l. 

4.2. The best model of the surface brightness profile 

We fit the model emissivity profile from the Bondi accretion 
flow convolved with instrument response and with the PSE to 
the observed surface brightness profile. In Eig. |3] (lower panel) 
the residuals of fitting up to 10" are shown. The fit is not sta¬ 
tistically acceptable, showing that the uniform Bondi flow is not 
a good representation for the hot plas ma at distances up to 10" 
from S gr A*. In previous analysis bv IShcherbakov & BaganoS 
(1201 Ol) a dynamical model of the hot plasma interacting with 
stars was fitted to the surface brightne ss profile up to the 5" from 
GC. Additionally. IWang et al.l(l2013h pointed out that Bondi ac¬ 
cretion flow ends at the radius ~ 3.6 x lO^Rschw which is 3.8" in 
the case of M - 4.4 x 10 ®Mq black hole. 

Eor this reason, we restricted our modelling to smaller radii 
and obtained a good fit up to 3" from Sgr A*. The best theoret¬ 
ical fit profile is shown with a dashed line in the upper panel of 
Eig.|6] and the parameters of the spherical accretion model are: 
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Table 3. List of spherical accretion model parameters (ng oo and Tea) with the corresponding mass accretion rate (Eql2]l, the speed of 
sound at the infinity, and Faraday Rotation measures for assumed radial, coherent, moderate (J3 - 1) and weak (J3 - 10^) magnetic 
fields. Fifth and sixth columns show RMgc computed towards Galactic Center, while seventh and eight columns represent RMpuisar 
computed towards pulsar PSR J1745-290 0, located 3" away fr om Sgr A* in the projected distance. For co mparison, the measur ed 
value of RM^^ = 5.6 + 0.7 x 10^ rad/m^ (iMarrone et al.llTOOTh . and = 6.696 + 0.004 x 10"^ rad/m^ (lEatough et al.ll2013h . 


Poo 

Too 

M 

Cs,oo 

RMgc 0? = 1) 

RMgc{jS= 10') 

RMpulsar {P - V) 

RMpa^ar iP = 100) 

[cm“^] 

[keV] 

[Moyr-‘] 

[ems-*] 

[rad/m'] 

[rad/m'] 

[rad/m'] 

[rad/m'] 

1 

1 

4.6 X 10-' 

3.9 X 10' 

4.4 X 10** 

1.4 X 10® 

4.5 X 10' 

40 

1 

3.5 

7.1 X 10-** 

7.4 X 10^ 

2.6 X 10' 

8.4 X 10' 

6.6 X 10' 

70 

1 

4 

5.8 X 10-** 

7.9 X 10'' 

2.0 X 10' 

6.2 X 10' 

7.1 X 10' 

70 

1 

5 

4.1 X 10-** 

8.9 X 10' 

1.2 X 10' 

3.8 X 10' 

7.8 X 10' 

80 

1 

6 

3.1 X 10-** 

9.7 X 10'' 

7.9 X 10® 

2.5 X 10' 

8.4 X 10' 

80 

1 

8 

2.0 X 10-** 

1.1 X 10** 

4.1 X 10® 

1.3 X 10' 

9.6 X 10' 

1.0 X 10' 

1 

16 

7.2 X 10-** 

1.5 X 10** 

9.0 X 10® 

3.0 X 10' 

1.3 X 10' 

1.3 X 10' 

18.3 

3.5 

1.3 X 10-" 

7.4 X 10' 

2.1 X 10'' 

6.6 X 10® 

5.2 X 10^ 

5.2 X 10' 

18.4 

3.5 

1.3 X 10-^ 

7.4 X 10' 

2.1 X 10" 

6.6 X 10® 

5.3 X 10" 

5.3 X 10' 

35 

1 

1.6 X 10-*’ 

3.9 X 10' 

9.1 X 10*" 

2.9 X 10' 

9.3 X 10" 

9.3 X 10' 

35 

3.5 

2.4 X 10-** 

7.4 X 10' 

5.5 X 10" 

1.7 X 10® 

1.4 X 10® 

1.4 X 10" 

35 

4 

2.0 X 10-® 

7.9 X 10' 

4.1 X 10" 

1.3 X 10® 

1.5 X 10® 

1.5 X 10" 

35 

5 

1.4 X 10-** 

8.9 X 10' 

2.5 X 10" 

7.8 X 10® 

1.6 X 10® 

1.6 X 10" 

35 

6 

1.1 X lO-*" 

9.7 X 10' 

1.6 X 10" 

5.2 X 10® 

1.7 X 10® 

1.7 X 10" 

35 

8 

7.2 X 10-’ 

1.1 X 10** 

8.6 X 10** 

2.7 X 10® 

2.0 X 10® 

2.0 X 10" 

35 

16 

2.5 X 10-’ 

1.5 X 10** 

1.8 X 10** 

5.8 X lO'* 

2.8 X 10® 

2.8 X 10" 

60 

1 

2.8 X 10-” 

3.9 X 10' 

2.0 X 10" 

6.4 X 10' 

2.1 X 10® 

2.1 X 10" 

60 

3.5 

4.2 X 10-® 

7.4 X 10' 

1.2 X 10‘® 

3.9 X 10® 

3.1 X 10® 

3.1 X 10" 

60 

4 

3.4 X 10-** 

7.9 X 10' 

9.1 X 10" 

2.9 X 10® 

3.3 X 10® 

3.3 X 10" 

60 

5 

2.5 X 10-® 

8.9 X 10' 

5.5 X 10" 

1.7 X 10® 

3.6 X 10® 

3.6 X 10" 

60 

6 

1.9 X 10-® 

9.7 X 10' 

3.7 X 10" 

1.2 X 10® 

3.9 X 10® 

3.9 X 10" 

60 

8 

1.2 X 10-® 

1.1 X 10** 

1.9 X 10" 

6.1 X 10® 

4.4 X 10® 

4.5 X 10" 

60 

16 

4.6 X 10-’ 

1.5 X 10** 

4.1 X 10** 

1.3 X 10® 

6.2 X 10® 

6.2 X 10" 



Fig. 6. Surface brightness profile up to 3" away from Sgr A*. 
Red points represent observed profile, and are connected by con¬ 
tinuous line. The best fit model with T™' = 3.5 + 0.3 keV and 
= 18.3 + 0.1 cm“^ is presented by blue short dashed line, 
and the fit residuals are presented in the lower panel. 


= 18.3+0.1 cm“^ and T™' = 3.5 ±0.3 keV. The correspond¬ 
ing integrated mass accretion rate is M - 1.3 x 10“® Mq yr“' 

. The reduced of our best fit is 1.19. Residuals are presented 
in the lower panel of Fig. All data points are within 2 sigma 
from the model. 


From the best fitting model we can estimate profiles of tem¬ 
perature and density of the hot plasma. We present them in Fig. |7] 
with solid lines. As a result of our fitting we found the values of 
density and temperature at the distance of 1.5" to be We = 29.2 
cm“^ and Tg = 4.8 keV, re spectively. These values are dif- 
ferent than those obtained by IShcherbakov & BaganoftI d2010h 
(wg = 130 cm“^ and Tg - 2 keV). 

The mass accretion rate from our Bondi best fit model is 10- 
1000 times larger in comparison to the mass accretion rate limits 
set by the measurements of the Faraday rot ation near the cen¬ 
tral black hole 10~^ < M < 10“^ Mq yr ' (iBower et alj|2005l : 
IMarrone et al.ll206^ l2007t) . Nevertheless, the derived mass ac¬ 
cretion rate is always model dependent and the direct compari¬ 
son is meaningless. The above limits computedbv iMarrone et al.l 
(l2006l see Eq. 8 in this paper) strongly depend how the Faraday 
rotation measure (RM) relates to the accretion flow and magnetic 
field projected on the line of integration. To avoid any inaccura¬ 
cies, we present here RM computed from our models according 
to the formula; RM=8.1 x 10^ J rigEdl, where electron density rig 
is in cm~^, the path length dl is in pc, and the magnetic field B 
- in Gauss (iGardner & Whiteoakll 19661) . the last defined by stan¬ 
dard formula B^ = 87rPgas//i, wherey6 is the plasma parameter. 

There is no direct constraint on the magnetic field value in 
Sgr A*, therefore we consider the case of moderate, and very 
weak magnetic field. We present RM computed towards the GC 
from all our models in Tab. [3] Fifth column is for P - i (mag¬ 
netic field in equipartition with the gas), while sixth column 
is for p - 10^. For p — 1, RMgc for the best fitted model 
does not agree with measured value RM^^, = 5.6 ± 0.7 x 10^ 
rad/m^ (IMarrone et alJlTOO^ . But it fully agrees, if we assume 
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Fig. 7. Temperature and density profiles around Sgr A* for our 
best fitted model T™' = 3.5 + 0.3 keV and = 18.3 + 0.1 

—T, ^ ^ 

cm“^, upper and lower panel respectively. The asymptotic val¬ 
ues of 7’°“'^ and are marked as a dashed lines. Additionally, 
we show the radial dependencies of the temperature and density 
in RIAF model without (s=0) and with (s=l) an outflow. The 
normalization of RIAF models is arbitrary. 


that the magnetic field in the radial direction is very weak i.e. 
P - 10^. We therefore admit, that our best fit model is not well 
constrained within approximately lOORschw or the magnetic field 
in radial direction is very weak. When we start the integration of 
our model from 5 x 10“^" (which corresponds to the distance of 
464Rschw), the resulting RMgc = 6.91 x 10^ rad/m^ with mag¬ 
netic field being in equipartition. In this case, our result is in quit 
go od agreement with ob servations, especially that in the paper 
by iMarrone et al.l (l2006h (Fig. 4) the integration of models was 
typically done from 300 Rschw, and than comparred to the data. 

For parameter p - 10^, the magnetic field in the center of 
the flow is B = 5 X 10“^ Gauss, which is not in agreement 
with recent found that the streng th of B in GC is around hun¬ 
dred Gauss (lEatough et alj|2013l) . Nevertheless, this value was 
given by authors after they have measured RM toward the pul¬ 
sar PSR J1745-2900 near the GC. Since all model computations 
are not precise very close to the black hole, we decided to com¬ 
pare the modelled RM towards the above pulsar, located 3" in 
the projected distance from the black hole. Such value was re- 
ported recently to be equal RMp^j„^ = 6.696+0.004x10"^ rad/m^ 
dEatough et al.ll2()T^ . We project our density profile and mag¬ 


netic field on the direction towards the pulsar and compute mod¬ 
elled RMpuisar in case of p - \ and p - 100. We present results 
in Tab. |3] seventh and eight column respectively. 

The RMpuisar for the best fitted model agrees very well with 
the observed value for p - Eor such value of p parameter 
the magnetic field towards Sgr A* is B = 157 Gauss, from our 
best fitted model. We conclude, that our model density and tem¬ 
perature profiles, and hence mass accretion rate of about 10 ® 
Mq yr ' around 3" are consistent with RM measured toward the 
GC pulsar. It is possible that the mass accretion rate is reduced at 
distances closer to the black hole. However, these outflows can¬ 
not be accounted for in our simple spherically symmetric model 
that is fitted to spherically averaged brightness profile. 

Eor comparison Eig.Qalso shows the radial profiles of tem¬ 
perature and gas number density from a nalytical radiativelly 
inefficient accre tion flow (RIAE) models (IShaniro et alJll97ra 
iRees et al.lll982h . In such models the density profile is given by 
n ~ where s is parameter that if positive then the RIAE 

launches an outflow. In Eig.|6]lower panel, we plot density profile 
for s - 0 (no outflow) and s - 1, by big and small dots respec¬ 
tively. The density normalization is arbitrary. In RIAE the tem¬ 
perature profile follows T ~ r'® where parameter 0 < 0 < 1 (the 
same figure upper pane l). Typically in RI AE 6 - \ , i.e. tempera¬ 
ture has a virial profile (IWang et alj2013h . In Bondi spherical ac¬ 
cretion below the sonic point the density profile follows the same 
radial dependency as in RIAE without the outflow (i.e. s=0), and 
the temperature is virial. However, the Bondi solutions at large 
radii naturally connects to constant values, and both, density and 
temperature profiles deviate from the n ~ and T ~ 
laws, despite of the fact that there is no outflow in this solution. 

Our result does not contradict other accretion flow models 
considered in case of Sgr A*. It shows the position of stagnation 
radius to be at 3". Outside this distance, profiles derived by us 
may not work. 

5. Conclusions 

In this paper we presented 134 ks Chandra ACIS-I observa¬ 
tions of the GC containing the Sgr A* region. After analysing 
flux images in the continuum (0.5 - 8 keV) and iron line (6.3 
- 7.0 keV) bandpasses, we have extracted spectrum of 5" cir¬ 
cular region around Sgr A*. Our spectral model consists of a 
thermal plasma for the continuum plus Gaussian profiles for the 
Kq, emission of Ee, S, Ar and Ca lines. The continuum is well 
fitted by a thermal bremsstrahlung emission indicating tempera¬ 
tures for the plasma around Sgr A* in the 2.2 - 2.7 keV range, 
depending on the choice of background. The corresponding ab¬ 
sorbing column densities towards the source are 10.5 - 9.2x10^^ 
cm“^, and the corresponding EWs of Ee^:^ line are 1.19 - 0.91 
keV. Our spectral fitting indi cates lower t emper ature by about 
0.5 keV than this obtained bv lWang et al.l (1201 3l) with the same 
model. Nev ertheless, the streng th of iron line consistent with that 
reported bv IWang et al.l (l20 13 ) at the Icr confidence level. 

Eurthermore, we studied the brightness profile of the hot 
plasma in the vicinity of the Sgr A* black hole. We found that 
the classical Bondi accretion can reproduce observed brightness 
profile up to 3". This result implies strong constrains on the po¬ 
sition of stagnation point always considered to be present in the 
complicated flow around Sgr A*. The best fit model for surface 
brightness profile requires that the hot plasma outside the flow to 
have temperature = 3.5 +0.3 keV and density n^ = 18.3 +0.1 

-3 

cm . 

Sgr A* has been r ecently resolved down to 1.5" and studied 
by I Wang et ^ d2013l) . They have reported an X-ray image of 


9 















































A. Rozanska R Mroz M. Moscibrodzka M. Sobolewska T. P. Adhikari: X-ray observations of the hot phase in Sgr A* 


the source based on a 3 Msec Chandra observation. The best bt- 
ted model was RIAF (radiatively inefficient accretion flow), but 
they did not show any resulting temperature and density profiles. 
Only a general description for the temperature dependence on 
the distance was given as: T oc where 0 > 0.6. Additionally, 
they concluded from the weak line of H-like iron Fe XXVI, that 
the amount of gas at temperatures ~ 9 keV is very small. Our 
temperature profile predicts the temperature of 9 keV at the dis¬ 
tance 0.4" from the black hole. Such plasma would have number 
density equal to 70 cm“^. Both temperature profiles, RIAF and 
Bondi, agree up to 0.1", as it is seen in Fig. [T] But Bondi flow 
fitted by us predicts higher density by one order of magnitude. 

All models show different temperatures at the outer parts 
of the flow Tg., = 1, 2, and 3 .5 keV respectively in the 
RIAF model (IWang et al. 20j^, hot plasma fed by stars 
dShcherbakov & BaganoffiboiOll . and the only Bondi flow re- 
port ed in this paper. A dditionally, spectral fitting done by us 
and IWang et ^ (1201 3h indicates slightly different temperatures 
(2.7 and 3.5 respectively). Those results show a big complex¬ 
ity of matter structure around Sgr A*, and additional indica¬ 
tions should be taken into account when study ing hot plasma 
in this region. For instance IWang et ^ (l2013l) also found the 
He-like Fe Kq, line. Therefore, we argue, that up to 3" Bondi 
flow works, but further away hot matter may be in the form 
of two-phase medium. We speculate that the Fe Kq. line emis¬ 
sion at 6.69 keV can be attributed to the reflection from clouds 
of the warm gas. This illuminated radiation originates from 
the very innermost nucleus and depending on the luminosity 
state of Sgr A*, can provide strong dumpin ess and trigger 
forma tion of such hot clouds, as proposed by iRozanska et al.l 
(|2014|) . However, it is critically important to confirm that the 
Bondi model fitt ed here does not o verpredict the strength of 
the FeXXVI line. IWang et ^ (1201 3h concluded that due to the 
weakness of FeXXVI line, there must be a little gas at tempera¬ 
tures larger than 9 keV. In our Bondi model the gas temperature 
reaches 9 keV at 0.4” (compared to 0.1” in the RIAF model). 
If the strength of this line computed from our best fitted model 
of Bondi flow disagrees with the observed one, this may rule out 
of our model despite its agreement with the surface brightness 
profile. To fully check our predictions we plan to construct the 
model of thermal plasma. Also, better spectrum will be needed 
for this kind of analysis. 

The mass accretion rate from our Bondi best fit model is 
1.3 X 10“® Mq yr“'. This is out of the range of accretion 
rates predicted from m easurements of the Faraday rotation (RM) 
(iMarrone et al.l 1200^ . Nevertheless, we point out that the de¬ 
rived mass accretion rate limits are always model dependent. 
Especially we can modify them by changing the structure and 
strength of a magnetic field. For our best fitted model the com¬ 
puted RMgc agrees with observed one toward Galactic Center 
with assumption of very weak radial magnetic field. Moreover, 
RMpuisar computed from the best fit model towards the pulsar 
PSR J1745-2900 located at 3" projected distance from GC is 
in very good agreement with observations assuming moderate 
magnetic field, /3 - 1. This result clearly shows that any limi¬ 
tation on accretion rate from Faraday rotation strongly depends 
on the model and orientation of magnetic field. We admit that 
the best way to compare the accretion model to Faraday rotation 
measure is not by comparing accretion rates, but by self consis¬ 
tent calculations of RM from the model. 

Profiles of temperature and density resulting from our fitting 
put constraints on the physical conditions in the vicinity of the 
GC. It may provide information on the value of the gas pressure 
of the hot plasma. Together with radiation pressure taken from 


observed luminosity states of Sgr A* we can estimate the size of 
the t hermal instability and eventual cloud formation in this re¬ 
gion (IRozanska et al.l2014l) . Such mul ti-phase co mplex structure 
is observed in our GC (for review, see lCoxll2005h . This structure 
is important from the point of view of feeding the central black 
hole as well as the possibility of an accompanying outflow. 

To study any possible outflow related tot he Sgr A* accre¬ 
tion, we should search for spatially resolved regions of the hot 
plasma farther away from the Sgr A* region, where the line 
emission is prominent. For example, strong iron Kq. line emis¬ 
sion detected in our observations occurs in the Sgr A East re¬ 
gion. The origin of this emission is strongly affected by SNR, 
and more advance spectral model should be used to describe the 
emerging spectra. We address these studies in the forthcoming 
paper. 
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